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The method of integral transforms is first applied for studying the '^He 
longitudinal response functions. The transforms are calculated from localized 
bound-state-type solutions to an inhomogenous Schrodinger-type three-body 
equation. The multipole expansion of the solutions is considered. The whole 
calculation is checked using a sum-rule test. Several versions of local s-wave 
spin-dependent potentials supplemented with a singlet p-wave potential and 
with the proton-proton Coulomb interaction are used as a two-nucleon in- 
put. The conventional charge density operator involving free-nucleon form 
factors is utilized. The three-body equations are solved with the help of the 
hyperspherical expansion and a complete convergence is achieved. The final- 
state interaction is thus fully taken into account. The q = 300 MeV/c and 
500 MeV/c values are considered. It is found that the contribution of the 
T = 3/2 final states to the problem is suppressed and it amounts about 15%. 
This might be ascribed to symmetry requirements. The contributions of the 
p-wave NN interaction and of the Coulomb interaction are found to amount 
several per cent. Uncertainty due to different choices of s-wave NN forces is 
of a similar magnitude provided that the low-energy A^A^ data are properly 
described. The results are compared with the integral transforms of the exper- 
imental response functions. For q = 300 MeV/c experimental and theoretical 
results coincide within their uncertainties. For q = 500 MeV/c a noticeable 
difference is detected. 
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Many calculations testify to the fact that the conventional form of the nuclear charge 
density is inapplicable for description of the elastic form factors of three- and four-nucleon 
nuclei at q > 2.5 fm~^ values. However the elastic scattering occurs with a quite low 
probability at such q values and some non-typical configurations including those where all 
the nucleons are close together may contribute substantially. In this connection it seems 
important to test a form of the nuclear 4-current in the inelastic processes and to study 
the (e, e') response functions. This requires a proper account of the nuclear final-state 
interaction. 

We study (e, e') response functions 

= Emo / < ^/l^l^o > \'^^^f - ^0 - e) (1) 
where e is the nuclear excitation energy. Direct evaluation of R with a full account of the 
final-state interaction is quite complicated even in the three-body case. Indeed, it requires 
obtaining the whole set of the final-state continuum wave functions ipf at the same energy 
and summing their contributions. As far as we know only one such an investigation was 
performed in the literature [1]. A local s-wave central potential [2] was utilized there while 
interaction in higher partial waves and Coulomb interaction were disregarded. 

In this paper we present the first microscopical study of the ^He longitudinal response 
functions Ri with the help of the method of integral transforms [3-6] . While the final-state 
interaction is fully taken into account in our approach we avoid calculating the continuum- 
spectrum wave functions ipf. We explore an ability of the method in an A=3 problem. We 
study sensitivity of the results to a choice of NN force and we clarify a role of an interaction 
in higher partial waves and that of the Coulomb interaction. We compare our results with 
experiment [7,8] at ? = 300 and 500 MeV/c. 

Define the reduced transition operator and the response function 

O = [GUQT'O, R{q,u) = [G%{Q)]-'R{q,u). 

Here = q^-u^ and [9] G'^(g) = [l + gV(4M2)]-V2G'^(g) where G% is the proton Sachs 
form factor. We calculate the integral transform of the response [3,10] 
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Hq,a)= {a + e)-'R{q,uj)de (2) 

instead of the response itself. We use the conventional single-nucleon expression for the 
charge density O. Then to a very good approximation one can disregard the cu-dependence 
of the O operator and use the expression 

where r'j = rj — Rem,. It has been shown [3,10] that a) can be calculated by first solving 
for the localized solution to the following inhomogenous equation 

{H~Eo + (t)§ = O^Po (4) 

where as in Eq. (1) i/jq is the ground-state wave function and Eq is the ground-state energy. 
In terms of \E' we have [3,10] 

^q,a)=<^\d^o>-cr-'Rei (5) 

where Rei is the elastic contribution to the response. 

The solution to Eq. is much easier to obtain than the functions ipf entering Eq. (|l]). 
Indeed, in contrast to the latter functions there is no need to impose the complicated large- 
distance boundary conditions in order to fix a solution. The only condition of vanishing of 
the solution at large distances suffices. Therefore methods that are used in solving bound- 
state problems can be utilized here. Below we use real a values in Eqs. (0), (^. 

We have two possible ways to connect our theoretical calculations with experiment. One 
way [10] is to compare a) with the same quantity obtained from the experimental R{q, e) 
using Eq. (2). Another way [3] is to consider Eq. (^) as the integral equation and invert 
to obtain theoretical R{q,u) and then compare the responses themselves. In the present 
work we use the first approach. This provides one only with a limited information since 
some R features are smeared out at performing the transform (|^). Nevertheless, this still 
allows testing the O transition operator in Eq. (|l|). If, for example, the nucleon form factors 
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should be modified inside the nucleus this would produce a change in $(g, a) of a magnitude 
comparable with the change in (G^)^. 

Concerning the inversion problem, another version of the method [6] better suits for this 
purpose. The many-body equation to be solved in this version is similar to Eq. (|^) yet. The 
only difference is that complex a values are utilized. It is shown below that the left-hand 
side of Eq. can be obtained from Eq. (1) with high accuracy and we achieved comparable 
accuracy for the complex a values as well. 

The calculations are performed at g = 300 and 500 MeV/c and they are compared with 
the Saclay [7] and Bates [8] experimental results. In this first calculation we use four versions 
of effective central local s-wave spin-dependent A^A^ potentials [2,11,12] supplemented with 
a realistic singlet p-wave NN potential [13] and with the proton-proton Coulomb interac- 
tion. Only the components of latter interaction which are diagonal in the isospin T =1/2, 
3/2 quantum numbers are retained in the calculation. Even these components prove to 
contribute little to the results, see below. 

Under these assumptions on the nuclear dynamics, Eq. (|) is split into independent sets 
of equations with a given orbital momentum L, and isospin T of the system. It is convenient 
to calculate the right-hand sides of these equations in the following way. Since ipo has L = 
then only the components ~ J2ml '^LAiLi^'j)'^LML(^) fro^i the expansion of exp(iqr^) from 
Eq. (^) contribute to the problem for a given L value. R is independent of a q direction 
due to averaging over Mq in Eq. (1). Let q be directed along the z axis. Then only the 
components with = give non-zero contributions and hence only the components of 
with Ml = are different from zero. We have 

iH-Eo + a)^LT = OltV-o, (6) 
d,T = E \TMt X TMr\ E[^^ + ^^^]j.(gr,)F.o(i-;), (7) 

^LT^q, ct) =< ^ltIOlt^o > -SLocr~^Rei, (8) 

oo 

$(g,a)=47rE(2^ + l) E '^LTiq,cj). (9) 

L=0 T=l/2,3/2 
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The functions \1/lt have the same spin = 1/2 as ipo- 

We solve Eqs. (H) using the hyperspherical expansion. Eqs. (6) turn into sets of algebraic 
equations. Denote K the hyperspherical momentum and [/] the type of symmetry of the 
spatial components of the basis functions. Our basis functions are of the form 

Rn{p)^^KLMj;^SMsTMT]i(^^ Cr^j, T^j) . (lO) 

Here p is the hyperradial variable and fl are hyperangular variables. The index i enumerates 
the functions with the same other quantum numbers. We use the Laguerre-type hyperradial 
basis functions [14] . We have developed a computer code to construct complete sets of the F 
basis functions from Eq. (10) with arbitrary quantum numbers using the the Raynal-Revai 
transformation [15]. The coefficients of this transformation are evaluated using the recurrent 
formula of the K ^ K + 2 type [16]. 

The L values up to Lmax = 8 are retained which provides convergence in Eq. (9). The 
net number of the basis functions with the same K, L, S, T and [/] values grows linearly 
with K which in general leads to systems of linear equations of high dimensions. But in 
fact these dimensions can be considerably reduced. For this purpose we use the fact that 
nuclear forces we utilize here only act in the two A^A^ orbital states. Namely, we specify 
orthonormalized basis functions in such a way [17] that only for two functions from each set 
of functions with the same K, L, S, T and [/] the matrix elements of A^A^ force are different 
from zero. Only such functions are coupled together in Eqs. (6). Other functions do not 
include s or p components for each nucleon pair. They lead to decoupled equations in Eqs. 
(6) and they correspond to free-motion final states. 

We retained the basis functions (10) with K < Kmax = 30 in our calculations. This 
provides practical convergence of the results. The convergence trends are illustrated in Fig. 
1 where dependence on K^ax of the calculated $l=2,t=i/2(q'5 cr) contribution is shown at 
q = 300 MeV/c and a = 1 MeV as an example. In order to check stability of our results 
with respect to K^ax we performed their extrapolations to K^ax = oo. This was done using 
three-parameter fits of the form ^{Kmax) = ^{oo) — c/K'^^^. Extrapolated values $(oo) 
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differ from ^{Kmax = 30) less than by 1%. Convergence with respect to hyperradial basis 
functions from Eq. (0) is achieved, too. The maximum number of coupled basis states 
equaled 861 at solving Eqs. (6). The calculated transforms are smooth functions similar to 
shown in Fig. la from Ref. 4 for the deutron case. 

We make a comment concerning calculation of $ at small a values. It is necessary to 
avoid large cancellations in the right-hand side of Eq. (H). This is achieved if one uses in 
Eq. for L = 0, T = 1/2 the same Kmax value as that at calculating ipQ. Then the pole 
terms cancel exactly. 

There exists a test which enables us to check the calculation as a whole. Namely, the 
leading term of a) at high a values behaves as a~^. This term can be calculated 
independently using the sum rule, 

/oo „ — „ _ „ 

R{q,u;)de = ^^^^ < iJo\0^0\^o > -Rei- (H) 
-min 

This allows one to check the right-hand side of Eq. (|^). Besides checking the calculation the 
test allows one to verify whether at high a values the results are stable against increasing 

KfYiax &nd -^^max-Q 

Various contributions to Eq. (9) are just the integral transforms of the corresponding 
contributions to R. We note that the T = 3/2 contributions prove to be suppressed in 
comparison with the T = 1/2 ones. The net relative T = 3/2 contribution to a) ranges 
between 12 and 17% for all the a values and the two q values considered. The contributions 
to the initial responses R should be of a similar magnitude. The reason for the suppression 
may be in that the spatial component of the final-state wave function which are symmetrical 



This test, however, is not applicable for checking the matrix elements of the Hamiltonian in Eqs. 
(6). (When a tends to infinity these matrix elements can be neglected.) In a preliminary version 
of this work [18] a considerable disagreement with experiment was found. This was due to an 
algebraic error. For mixed symmetry states the matrix elements of kinetic energy were erroneously 
twice as large as correct ones. 
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with respect to particle interchanges are present at T = 1/2 only. These components may 
provide more internucleon attraction than those of the other symmetries which increases 
amplitudes of the final-state wave functions inside the reaction zone. 

We estimated an influence on the results of NN interactions in higher partial waves and 
of the proton Coulomb interaction. The results are shown in Fig. 2 for g = 300 MeV/c. We 
performed the calculation for the MT(I-I-III) potential supplemented with the realistic p- 
wave singlet interaction from Ref. 13 and with the Coulomb interaction. Then we switched 
off the p-wave interaction. Curve 1 shows the relative change in the results. The contribution 
from the triplet p-wave force is believed to be of the same size. After that we switched off the 
Coulomb interaction in addition. This produces Curve 2. It is seen that the contributions 
of the p-wave interaction and of the Coulomb interaction do not exceed several per cent or 
standard uncertainties in $ functions obtained from experimental data. For q = 500 MeV/c 
these contributions are even much smaller. 

Furthermore, we studied dependence of the results on the choice of the A^'A'" force. Central 
local s-wave A^A^ forces utilized in our calculation include the MT(I-I-III) [2], S2,S3 [11] 
and EH [12] potentials. The MT(I+III) and S2 potentials reproduce the A^A^ low-energy 
properties and s-wave A" A" phases up to high energies. The S3 potential fits the low-energy 
data and yields nearly correct values for the binding energies and rms radii of ^He and 
^He. The EH potential fits the s-wave A"A^ phases up to high energies but it does not 
reproduce properly the low-energy A^A^ data (see Ref. 11). As above all these forces were 
supplemented with the realistic p-wave singlet interaction [13] and with the proton Coulomb 
interaction in our calculations. Curve 3 in Fig. 2 represents the relative difference in $ 
between S2 and MT(I-I-III) potentials. In case of S3 and MT(I-I-III) potentials such a 
difference is approximately of the same value and of opposite sign. These differences do not 
exceed several per cent. Curve 4 represents such a difference between EH and MT(I-I-III) 
potentials. The latter difference exceeds ten per cent at small a values which corresponds 
to substantial differences in low-cncrgy parts of the responses. This can be attributed to 
the above-mentioned drawback of the EH potential. For q — 500 MeV/c differences in the 
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results for different potentials become considerably smaller. 

Before comparing our calculation with experiment we make some comments on obtaining 
the integral transforms, Eq. (0)), of the experimental responses. In order to perform the 
integration in Eq. (2) with a sufficient accuracy and in particular to estimate the contribution 
from the unavailable high-e tails of the spectra we approximate the experimental R functions 
by the following analytical expressions: a{uj—uJthreshY^'^ in the low cu region uJthresh < < ^i-, 
then Yl,n=o^n^^ in the region of the peak ui < u < U2 and Z]^=o Cn(^^2/t^)°'^" in the 
region beyond the peak. The parameters a, 6„, c„, a and uj2 are chosen from the least-square 
procedure with additional requirements imposed of continuity of the fitting spectra and 
their first derivatives at uji and 002 points. It turns out that a good description at quite wide 
ranges U2 < uj < uJmax of the spectra beyond the peaks are provided with a single term [19] 
~ (In case, say, exponentially decreasing tail-terms the description is worse.) The best 
a values range between about 4 and 5 for all the q values considered. Similar a values were 
found [20] in the "^He case. We extrapolate the fitted spectra beyond the oJmax values in 
order to take into account the contributions from the unavailable tails. In case of data from 
Ref. [7] these contributions proved to be quite small. They reach their maxima at high a 
values where they are between 1 and 2%. In case of q=500 MeV data from Ref. [8] we did 
not succeed in producing stable extrapolations. 

In Fig. 3 we compare the calculated $ values and the values deduced from experiment. 
The theoretical calculation was done MT(I+III) A^A^ potential [2] supplemented as above 
with p-wave singlet interaction and with the proton Coulomb interaction. The rela- 
tive differences between the theoretical and experimental $ values are shown. Taking into 
account 5% systematic uncertainties of the experimental data and also above- considered un- 
certainties of theoretical calculations one can say that that there is no significant difference 
between experiment and theory at g = 300 MeV/c. The deviation at q =500 MeV/c may 
be considered as significant. Detectable differences betweem experiment and theory at such 
q values for another set [8] of data were also obtained in Ref. [1] where the responses were 
calculated directly. They may be attributed to relativistic effects. 
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In conclusion, we applied the method of integral transforms for studying ^He longitudi- 
nal response functions. This requires solving for a localized solution to an inhomogenous 
Schrodinger-type three-body equation. We elaborated techniques for this purpose and we 
obtained accurate solutions using central local NN potentials and the conventional expres- 
sion for the nuclear charge density. We found that A^'A'" interaction in higher patial waves 
and the proton Coulomb interaction play a minor role in the problem. Uncertainties in the 
s-wave NN force proved to be not substantial as well provided that the low-energy A^A'"-data 
are properly described. We compared the calculated integral transforms with experiment 
and we found that at q =300 MeV/c they agree with each other within their uncertainties. 
For q —500 MeV/c noticeable deviations are found. The results obtained make it possible 
to solve the corresponding problem in the quite important a-particle case, in particular. In 
addition, we found that the final-state T = 3/2 contributions to the problem are suppressed. 

We are grateful to K. Dow and C. Marchand for providing us with the experimental data. 
Helpful discussions with G. Do Dang and Yu. E. Pokrovsky are gratefully acknowledged. 
The work was supported by the Russian Fund for Fundamental Research under Grant No. 
93-02-14405. 
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CAPTION TO FIGURES 

FIG. 1. Dependence of $Lr(g, cr) on Krr^ax at L = 2, T = 1/3, ? = 300 MeV/c and cr = 1 
MeV. 

FIG. 2. Relative differences in cr) between MT(I+III) potential [2] supplemented with 
p-wave singlet interaction plus Coulomb proton interaction and other NN inputs. Curve 1 
- p-wave interaction switched off. Curve 2 - Coulomb interaction switched off in addition. 
Curve 3 - a difference with the S2 NN force [11]. Curve 4 - a difference with EH NN force 
[12]. (In the last two cases p-wave singlet interaction and Coulomb interaction are included 
in the calculation.) 

Fig. 3. Comparison of the calculated transforms cr) with experiment. Curve 1 is for 

experimental data from Ref. [7] for q — 300 MeV/c. Curve 2 is for experimental data from 
Ref. [8] for q — 300 MeV/c. Curve 3 is for experimental data from Ref. [8] for q = 500 
MeV/c. 
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